Slow-Roll Suppression of Adiabatic Instabilities in Coupled Scalar Field-Dark Matter 

Models 
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We study the evolution of linear density perturbations in the context of interacting scalar field- 
dark matter cosmologies, where the presence of the coupling acts as a stabilization mechanism for 
the runaway behavior of the scalar self-interaction potential as in the case of the Chameleon model. 
We show that in the "adiabatic" background regime of the system the rise of unstable growing modes 
of the perturbations is suppressed by the slow-roll dynamics of the field. Furthermore the coupled 
system behaves as an inhomogeneous adiabatic fluid. In contrast instabilities may develop for large 
values of the coupling constant, or along non-adiabatic solutions, characterized by a period of high- 
frequency dumped oscillations of the scalar field. In the latter case the dynamical instabilities of the 
field fluctuations, which are typical of oscillatory scalar field regimes, are amplified and transmitted 
by the coupling to dark matter perturbations. 
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I. INTRODUCTION 

Cosmology has provided evidence of a dark physics sec- 
tor which is necessary to account for about 95% of the 
cosmic matter content Despite the success of the 
ACDM model to fit all cosmological observations, the 
existence of the dark energy phenomenon as well as its 
relation to the abundance and clustering of matter in the 
universe still pose puzzling questions. 

Models of interacting dark energy-dark matter have 
been proposed to address such problems. In this scenario 
dark energy is a fundamental scalar field which directly 
couples to matter particles. This allows for a dynami- 
cal solution of the so called "coincidence" problem, since 
independently of the initial conditions the scalar inter- 
action drives the dark energy-to-matter ratio toward a 
constant value (see e.g. @, H, Hi)- These models are 
inspired by string and supergravity theories, where the 
compactification of extra-dimensions in the low-energy 
gives rise to massless scalars coupled to matter fields 
with gravitational strength. Therefore a distinct fea- 
ture of this scenario is that matter particles experience 
a long-range scalar force and acquire a time dependent 
mass which cause violations of the Equivalence Princi- 
ple (EP). The tight bounds imposed by EP tests are 
usually avoided as consequence of other possible mecha- 
nisms. As an example Damour and Polyakov have shown 
that in String Theory the couplings between the dilaton 
and different matter fields can be dynamically suppressed 
Q. An interesting possibility has been proposed in the 
"Chameleon" model [7[ , where the mass of the scalar field 
is assumed to depend on the local matter density. In 
such a case fifth-force effects can be strongly suppressed 
on Solar System scales, thus avoiding EP bounds. An- 
other possibility has been explored in 0] where the au- 
thors consider a dilatonic field to be differently coupled 
to various matter species such that the system can natu- 
rally evolve toward a late time attractor solution where 



General Relativity is recovered. Non-minimally coupled 
models can successfully describe the background expan- 
sion of the universe as probed by supernova type la lu- 
minosity distance or the position of the Doppler peaks 
in the Cosmic Microwave Background anisotropy power 
spectrum (see e.g. [H, [l(|). However testing the forma- 
tion of structure in the universe more than standard cos- 
mological tests may provide a key insight on this class of 
models. In fact the scalar coupling contributes to mod- 
ifying the clustering properties of matter, implying that 
an accurate study of the evolution of density fluctuations 
both in the linear and non-linear phase of collapse can 
identifying unique signatures of dark sector interactions 
PH - In the context of linear perturbation theory sev- 
eral interacting scalar field-dark matter models have been 
studied in the literature (see e.g. [3). In some specific 
realizations it was found that the growth of linear den- 
sity perturbations is spoiled by the presence of danger- 
ous instabilities [H, [Tj], as in the case of "Mass Vary- 
ing Neutrino" (MaVaN) models jla |. Recently a num- 
ber of works have analysed the stability of perturbations 
in more general setups. For instance in [161 ] the authors 
have studied models with a background evolution charac- 
terized by an adiabatic regime, and shown that unstable 
growing modes of the perturbations exist for couplings 
much greater than gra vitational strength. On the other 
hand the authors of [17l.[l8| have considered the case of an 
interacting dark energy component with constant equa- 
tion of state and found that for couplings proportional to 
the dark matter density the perturbations are unstable. 

In this paper we provide a more detailed study of these 
instabilities, particularly in relation to the specificities of 
the background scalar field evolution. The paper is orga- 
nized as follows: in Section [IT] we introduce the interact- 
ing scalar field-matter model as well as the background 
and perturbation equations; in Section [TTT1 we present the 
results of our analysis; finally in Section IIVI we present 
our conclusions. 
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II. INTERACTING SCALAR FIELD-DARK 
MATTER MODEL 

Let us consider a scalar field (j> with direct coupling to 
matter particles via a Yukawa term f(<f>/Mpi)tpip, where 
/ is the coupling function and ip is a Dirac spinor rep- 
resenting the matter field (Mpi = 1/ \/87rG is the re- 
duced Planck mass with G being the Newton constant). 
The effect of the scalar-dependent coupling is to induce 
a time-varying mass of the matter particles, hence caus- 
ing a violation of the EP. As mentioned in the previous 
Section, there are several ways to evade the tight bounds 
from EP tests. Here we assume that the scalar field only 
couples to dark matter particles. Therefore for the pur- 
poses of our analysis we neglect the baryon contribution 
and focus only on the cosmological evolution of the cou- 
pled scalar field-dark matter system. 

As in the case of the Chameleon cosmology [l9j], we 
assume the 0-field to have a self-interaction potential of 
runaway type in the form of inverse power-law: 
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where M is a mass scale and a is a positive constant. 
We consider a coupling function of dilatonic type, f(<p) = 
exp (/3(j)/Mpi), with (3 a dimensionless coupling constant. 
The background evolution of this system has been stud- 
ied in detail in @. 



FIG. 1: Scalar field effective potential at z = 10 3 ,10, 3 and 
(solid lines) for a — 0.2 and /3 = 1. The amplitude of 
the scalar potential M is set such that today Qdm = 0.24 
(yi<f, = 1 — Q.dm)- The dashed line corresponds to the position 
of the minimum of the effective potential at different epochs. 



A. Background and Linear Perturbation Equations 

Let assume a flat Friedmann-Lemaitre-Robertson- 
Walker metric (ds 2 — ~dt 2 + a(t) 2 dx 2 ), the evolution 
of the scalar factor is given by: 



pDM + tf/Z + V{(j>) 



(2) 



where pdm is the dark matter density and we have 
adopted Planck units (Mpi = I). The total energy 

momentum tensor of the system is conserved, Tufff = 



rpll(DM) rpH 



In contrast the non-minimal cou- 



pling implies that the energy momentum tensor of each 
individual component is not conserved. In such a case 
we can consider 



rpfl(DM) 

from which we obtain: 



rp^(DM) 



PDM + SHp DM = P4>PDM, 

4> + 3H<f>+V t<j , = ~(3p DM - 



(3) 
(4) 

(5) 
(G) 



Without loss of generality we can rescale the coupling 
function f(cj>) to its present value, f(<fio). Hence the so- 



lution to Eq. ([5]) is 



n<5 



(7) 



where p^ M is the present matter density. We may notice 
that as consequence of the scalar interaction the dark 
matter density deviates from the standard scaling a~ 3 . 
Furthermore for coupling values (3 > 0, the system of 
Eqs. (O-© describes an energy transfer from the 0-field 
to dark matter. In such a case the scalar field evolves in 
an effective potential 



KffW = V(4>) 



(0) 

PDM J3(4>-4>a) 



(8) 



which is characterized by the presence of a minimum. 

In figure Q] we plot the effective potential for (3 — I and 
a = 0.2 at redshift z = 1000, 10, 3 and respectively. 
For this choice of the model parameters we have (f>o ~ 
0.7605 as obtained by integrating numerically the system 
of Eqs. (H])-©. The dashed line in Fig. [1] corresponds to 
the position of the minimum at different epochs. 

In synchronous gauge the linearized equation for the 
dark matter density contrast 8dm, velocity gradient 
&dm, and field fluctuation S<fi are given by: 
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f38<t>, 
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8tj> + 3H5cj) 

= —PpdmSdm, 
where h is the metric perturbation given by 
2k 2 rj 8ttG 



h = 
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and 



V = 



PdmOdm + ak' 



(10) 



(11) 



(12) 
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In Section IIIII we will present the results of the numer- 
ical integration of this system of equations. However for 
a qualitative understanding of the conditions which lead 
to the onset of instabilities during the growth of the den- 
sity perturbations, it is useful to introduce an effective 
unified fluid description. 



B. Effective Unified Fluid Description 

The conservation of the total energy momentum ten- 
sor allows us to describe the interacting scalar field-dark 
matter system as a single unified fluid. The equation for 
the background density is given by 



Pt 



-3H(1 + w t )pt, 



(14) 



with pt — <p 2 /2 + V(4>) + pdm and wt = Pt I Pt, where 
Pt = <A 2 /2 + V ((/)). Similarly at linear order the pertur- 
bation equations in synchronous gauge read as: 

S T = —3H(c sT 




= -H(l 



where c 2 T = pt I Pt is the square of the adiabatic sound 
speed of the unified fluid and c 2 T = Spt / SpT is the square 
of the speed at which pressure perturbations propagate 
in the fluid rest frame. For a barotropic fluid with a 
constant equation of state (e.g. matter, radiation) c 2 s = 
c 2 = w. This is not the case for a generic fluid (e.g. 
scalar field), for this reason we may expect the effective 
unified fluid to be non-barotropic, (i.e. c 2 T ^ c 2 T ^ wt)- 
In terms of the scalar field and dark matter variables we 
have 



3H<P 2 



\2V d 



(3p 



DM \ 



3H<j> 2 + ZHp D M 



V,pS(j) + pdmSdm 



(17) 



(18) 



These relations provide us with a simple way of deter- 
mining the properties of the perturbation in the coupled 
system. For example in a given background regime insta- 
bilities of the perturbations may develop if the adiabatic 
sound speed acquire sufficiently negative values. 



III. SCALAR FIELD DYNAMICS AND 
EVOLUTION OF DENSITY PERTURBATIONS 

The non-minimally coupled scalar field model de- 
scribed in Section |TT] is characterized by the existence of 
an attractor solution which is set by the minimum of the 
effective potential. The minimum is given by V'g = 0, 
thus along the attractor solution the following condition 
is always satisfied: 



Va 



(0) 

PDM r 0(4>-<t>o) 

n3 ^ 



(19) 



Evaluating the derivative of Eq. {T]) and substituting in 
Eq. (TIT))) we obtain the time evolution of the field at the 



minimum: 



a+l 



= /3(0mi„-0 O ) 



(20) 



which depends on both the slope a and the coupling (3. 
Equation f|20[) is a non-linear algebraic equation which 
can be solved numerically through standard bisection 
methods (see dashed line in Fig. [1]). 

The field may reach the minimum from two different 
sets of initial conditions: <fii n i < 4>mi n (small field) or 
<\>ini > <^min (large field). In the former case <fi evolves 
over the inverse power-law part of the effective potential, 
where it minimizes the potential by slow-rolling as shown 
in Q. In fact one can easily verify that throughout the 
cosmological evolution the field mass (to 2 = ) as 
well as the ratio of its kinetic-to-potential energy satisfy 
the conditions to > H and <j> 2 /2V < 1 respectively. In 
contrast starting from large field values, <f> rolls towards 
the minimum along the steep exponential part of V^ff (0). 
Thus it rapidly acquires kinetic energy which subse- 
quently dissipates through large high-frequency damped 
oscillations around the minimum. 

As we shall see next, the growth of linear perturbations 
in these two regimes is significantly different. 



A. Adiabatic Regime 



As mentioned in Section III Bl we can obtain a qualita- 
tive insight on the stability of the perturbations in the 
coupled system by considering the effective unified fluid 
description. Let us evaluate the adiabatic sound speed 
Eq. (fl~7|) along the adiabatic solution Eq. (|19p : after ne- 
glecting the term proportional to the kinetic energy of 
the scalar field we have 



r 2 



-P 



3H' 



(21) 
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since <j> > it then follows that c 2 T < 0, implying that 
adiabatic instabilities may indeed develop. However we 
should remark that during the adiabatic regime the field 
is slow- rolling (i.e. 3Hcf> w 0), hence the term <f>/3H can 
be negligibly small compared to (3, such that c 2 T w 0, 
hence leading to a stable growth of the perturbations. In 
contrast instabilities will occur if the coupling assumes 
extremely large values, f3 ^> 3H/<j>. This is consistent 
with the analysis presented in [lj|, where the authors 
have suggested that during the adiabatic regime pertur- 
bations suffer of instabilities provided that f3 3> 1. Here 
we want to stress two main points which were not ad- 
dressed in that study: first of all that the rise of instabil- 
ities is suppressed by the slow-rolling of the field in the 
adiabatic regime, and secondly that exactly because of 
the slow-roll condition, instabilities can spoil the growth 
of dark matter perturbations only for large unnatural 
values of the coupling. To give an example let us as- 
sume that for a given model along the adiabatic solution 
the following condition occurs: c/)/3H <~ 10~ 2 . In such 
a case instabilities will develop only if the coupling con- 
stant P > 100, corresponding to a scalar fifth-force which 
is 2000 times greater than the gravitational strength. 1 
Moreover during the adiabatic evolution, Eq. HE 
reads as 



c sT 



1 J_ <>DM 

P &<t> 



(22) 



and assuming that the scalar field is nearly homogeneous, 
5(f> <C 5dm (in Planck units), we have c 2 T ~ /35<j)/SDM, 
and for f} ~ 0(1) this implies c 2 T « 0. In other words 
if the scalar field fluctuations are small with respect to 
the dark matter density contrast, then the coupled sys- 
tem behaves has a single adiabatic inhomogeneous fluid 



2 

sT 



0). 



These results are supported by the numerical study of 
the system of Eqs. (|9))- (IT3)) . with the scalar field evolution 
given by Eq. (j2TJ)) . We have set the model parameters to 
the the following values: a = 0.2, /3 = 1, with floM — 
0.24, H = 70 Kmr^pc -1 . As shown in [9] this model 
has the interesting feature that the background dynamics 
can mimic that of a phantom cosmology corresponding to 
an uncoupled dark energy model with slightly constant 
super- negative equation of state wjje = —1.1. 

The results of the numerical integration are shown in 
Figure [H In the upper left panel we plot the scalar field 
equation of state (solid line) and the equation of state 
for the effective unified fluid wt (dot line) . As we can see 
w<f, — —1, which is consistent with the fact that (j>/2>H 
is negligible, as can be seen from the plot in the right 



1 As consequence of the scalar interaction dark matter particles 
experience a gravitational force with effective Newtonian con- 
stant G e ff = G(l + 2/3 2 ). In contrast baryonic bodies may not 
experience such modification due to the non-linear nature of the 
scalar interaction 12011. 
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FIG. 2: Upper left panel: evolution of the scalar field equa- 
tion of state w$ and effective unified fluid equation of state 
wt', Right upper panel: evolution of the scalar field veloc- 
ity with respect to the Hubble rate; Lower left panel: redshift 
evolution of the adiabatic sound speed <? aT and propagation of 
pressure perturbations (? sT \ Right lower panel: Linear growth 
factor of the dark matter density contrast at k — 10 -3 , 10 -2 



and 0.1 Mpc" 1 



upper panel. We can also notice that the unified fluid at 
early times behaves as a matter component (wt = 0) and 
deviates toward negative values (—1 < wt < 0) as the 
0-field becomes the energetically dominant. In the lower 
left panel we plot the absolute value of c 2 aT , and c 2 sT {k) 
for three different scales k = 10~ 3 ,10~ 2 and 0.1 Mpc -1 
respectively. The adiabatic sound speed has negligible 
negative values and evolves with a trend that matches 
that of cb/3H, which is consistent with Eq. (f2"Tj) . We 
can also notice that the speed of propagation of pres- 
sure perturbations in the unified fluid remains is ~ 0. 
Hence during the adiabatic regime the interacting sys- 
tem behaves as a single inhomogeneous adiabatic fluid. 
In the lower right panel we plot the evolution of the dark 
matter density contrast normalized to the present value 
for k — 10~ 3 , 10~ 2 and 0.1 Mpc -1 respectively (for clar- 
ity we have displaced by a constant factor the different 
curves which would otherwise nearly overlap). As ex- 
pected these different modes manifest a standard power 
law growth and no instabilities are present. These re- 
sults have been obtained for an inverse power-law poten- 
tial, nevertheless they can be generalized to other scalar 
potentials, the only requirement is the existence of an 
adiabatic solution during which the slow-roll condition is 
satisfied. 
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B. Non-Adiabatic Regime: Large Field Oscillations 

Starting from initially large field values, the system 
evolves along a non-adiabatic solution characterized by 
rapid dumped field oscillations around the minimum of 
the effective potential. We can see this explicitely the 
upper left panel of Figure [3l where we plot the evolu- 
tion of the scalar field equation of state for the same 
model parameters as in Section IIII Al and obtained by 
numerically integrating Eq. © with initial conditions: 
<t>(a ini = 10- 5 ) = 0.15 > 0™ and fam = 0. We can 
infer the main features of the scalar field evolution from 
the behavior of it equation of state shown in the upper 
left panel of Figure [3] As we can see the field initially be- 
haves as a stiff component with = 1, this is because 
the field starts rolling on the steep exponential part of 
the potential, and consequently its energy is dominated 
by the kinetic term. As the field reaches the opposite side 
of the potential, it undergoes a series of high-frequency 
dumped oscillations around the minimum during which 
it dissipates most of its kinetic energy. It then sets on the 
inverse power-law part of the potential where it evolves 
along a tracker solution with ~ —2/ (2 + a) w —0.9. 

The evolution of density perturbations in the case of 
oscillating scalar fields has been widely studied in the 
literature, particularly in context of inflation 2j]]- From 
these studies it is well known that scalar field fluctuations 
are unstable during oscillatory regimes. In [22j the au- 
thors have presented a simple insightful explanation for 
the onset of such instabilities. The idea is to interpret 
the scalar fluctuation 8(f) as the separation between two 
particles whose dynamics is described by two coupled an- 
harmonic oscillators. Then a simple stability criterion is 
given by the relation between the frequency of the oscilla- 
tions ui, and their amplitude <j> (23j . Let us suppose that 
the frequency increases as the amplitude of the oscilla- 
tions diminishes, in such a case it has been shown that 
the distance between the two particles increases, thus 
causing the scalar field fluctuation to be unstable [22I ]. 
This is indeed what occurs in the interacting scalar field- 
dark matter system along the non-adiabatic solution we 
are considering. In fact we can see in the right upper 
panel of Fig. [3] that as the field starts oscillating, the 
frequency of the oscillations increases as the field ampli- 
tude diminishes, (du>/d(f> < 0). We can therefore expect 
the presence of instable modes. This is confirmed by the 
numerical solutions of 6<pk and 8dm obtained from the 
integration of Eqs. (|9"))- (fT5)) . The evolution of the scalar 
field fluctuation 5<pk is shown in the lower left panel of 
Fig. [3J We may notice the presence of an instability oc- 
curring roughly at the same time of the first oscillation, 
then followed by a second stage of exponential growth at 
the beginning of the second oscillation. From the plot in 
the lower right panel we can also see that the same insta- 
bility is passed to the dark matter perturbation, which is 
a direct consequence of the coupling terms in Eq. ^ and 
E q. (1101). Such unstable modes are similar to those found 
in [17lll8||. in fact by averaging over periods of time larger 




FIG. 3: Upper left panel: evolution of the scalar field equation 
of state w,),; Right upper panel: evolution of the scalar field; 
Lower left panel: evolution of the field fluctuations 5<j>h at k = 
10~ 3 , 10~ 2 and 0.1 Mpc -1 respectively; Right lower panel: 
evolution of dark matter density for fc-values as in the case of 
5(j>k- 

than the characteristic time of the oscillations, the scalar 
field behaves effectively as a dark energy component with 
a constant equation of state. 



IV. CONCLUSIONS 

We have studied the evolution of linear perturbations 
in the case of an interacting scalar field with runaway po- 
tential directly coupled to dark matter particles. We have 
specifically analyzed the stability of perturbations during 
the adiabatic evolution of the field, and shown that as 
consequence of the slow-roll condition the onsets of in- 
stabilities is largely suppressed. This can be explained in 
terms of the adiabatic sound speed of the effective uni- 
fied fluid. In fact during the adiabatic regime, despite 
being negative, it assumes negligibly small values and 
as consequence of this the growth of linear density per- 
turbations remains stable. On the other hand instabili- 
ties may develop in strongly coupled adiabatic regimes, 
with a coupling constant much greater than gravitational 
strength. Interestingly during the adiabatic evolution of 
the field the coupled system behaves as a single adiabatic 
inhomogeneous fluid. We have also shown that large in- 
stabilities can spoil the growth of linear perturbations in 
the case of non-adiabatic solutions characterized by large 
scalar field oscillations. It is well known that scalar field 
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fluctuations are unstable during oscillatory regimes, in 
such a case the scalar coupling amplifies and propagates 
such instabilities to the perturbations of the dark matter 
component. 

Our analysis suggest that under minimal natural model 
assumptions Chameleon-like cosmologies are not affected 
by instabilities of the perturbations and can provide a 
viable period of structure formation more than previously 
believed. 
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